Morphological and morphometric study of the scapulae of Korean wild deer

Korean water deer (Hydropotes inermis argyropus; Heude, 1884) and Siberian roe deer (Capreolus pygargus; Pallas, 1771) are Korean wild deer classified in the tribe Capreolini. C. pygargus in Korea were previously considered a single species; however, it was recently suggested that roe deer living on Jeju Island (Jeju roe deer; Capreolus pygargus jejuensis) is a distinct subspecies from roe deer living on the Korean peninsula (mainland roe deer; Capreolus pygargus tianschanicus) based on several studies demonstrating genetic and morphological features. In this study, we suggests that the scapular morphology and osteometric data can be used for interspecies discrmination between Korean wild deer. To compare the morphological characteristics of scapula among the three groups of deer, we analyzed the features and nine osteomorphological measurements of 31 H. i. argyropus (14 males and 17 females), 18 C. p. jejuensis (4 males and 14 females), and 23 C. p. tianschanicus (16 females and 7 males). The estimated ages of the deer were over 32–35 months. Data were analyzed by one-way repeated measures analysis of variance with post hoc Duncan test and discriminant functional analysis (DFA). H. i. argyropus and C. p. tianschanicus had the smallest and largest scapulae, respectively. The scapulae of the three Korean wild deer had a similar triangular shape, which was obscured by the tuber of the scapular spine, pointed acromion, broad infraspinous fossa, narrow supraspinous fossa, and partial ossification of scapular cartilage in older deer. H. i. argyropus had certain distinctive features, including a caudally pointed acromion, a notch between the supraglenoid tubercle and glenoid cavity (NBSG), a glenoid notch, and no sexual dimorphism, except for the longest dorsal length (Ld) and the scapular index (SI). C. p. jejuensis had a larger scapular index (SI) (61.74 ± 0.74%), compared with the SIs of H. i. argyropus and C. p. tianschanicus. The unique features of the scapula in C. p. jejuensis include its S-shaped cranial border. The C. p. jejuensis had a cranially pointed acromion, less frequent presence of glenoid notch and NBSG, short length of supraglenoid tubercle, and no sexual dimorphism. The C. p. tianschanicus had elevated cranial margin of the glenoid cavity, and frequent presence of glenoid notch and NBSG, similar to the H. i. argyropus. Similar to C. p. jejuensis, C. p. tianschanicus had a cranially pointed acromion. However, sexual dimorphism was observed in C. p. tianschanicus. DFA using osteometric data showed 97.22% of specimens were classified correctly into their species, meaning the osteometric parameters can be used for interspecies discrimination of Korean wild deer. Our findings indicate that the scapular morphologies of the three Korean wild deer have certain similarities and differences, suggesting that C. p. jejuensis are distinct from C. p. tianschanicus.

In this study, we analyzed the morphological features of scapulae in H. i. argyropus and roe deer. The animals used in this study were aged over 32-35 months. Left scapular morphology was analyzed using the methods described by Von den Driesch 25 .

Results
Scapular morphology. Scapular spine extended from the dorsal border to the neck of scapula. Dividing lateral side of the scapula into supraspinous and infraspinous fossae, the scapular spine was noted at a level near one-fourth of the length of its dorsal edge. The size of the spine gradually increased until the scapular spine tuber which was located at the proximal third of the length of the spinous process (Fig. 1). The free margin of the scapular spine was sinuously curved; caudally at the proximal third of the scapular spine, distally at the scapular spine tuber, and then cranially at distal third of it. The scapular spine tuber was visible as a slightly thickened free margin, and the proximal region of the tuber was weaker than its distal region (Fig. 1). The direction of the acromion differed between C. pygargus and H. i. argyropus. The acromion of H. i. argyropus exhibited a caudal curvature, except in a single male deer, whereas the acromion of C. pygargus was cranially curved (Table 1). Additionally, the acromion was sharp in all deer, and its tip ended at the level of the scapular neck (Fig. 2). Among the three Korean wild deer, H. i. argyropus had the sharpest acromion. The supraspinous fossa was smaller than the infraspinous fossa. Both fossae were triangular, narrowed toward the scapular neck, and merged at the acromion level (Fig. 1).
The cranial border of the scapula (i.e., the line connecting the cranial angle to the ventral angle) was sharper than the caudal border. Additionally, the free margin was curved toward the scapular neck. In H. i. argyropus, the cranial border was oriented cranially from the cranial angle, then oriented obliquely toward the scapular neck at the level of the scapular spine tuber (Fig. 1). After formation of the scapular notch, the cranial border is oriented cranioventrally toward the supraglenoid tubercle (Fig. 1). The cranial border of C. p. jejuensis scapula had a well-developed curvature, similar to a gentle "S" shape ( Fig. 1). It arose cranially from the cranial angle, was oriented along the scapular spine at an intermediate scapular height, and ended at the supraglenoid tubercle, making broad area of supraspinous fossa. In contrast, the cranial border of the Korean peninsula-dwelling C. p.  www.nature.com/scientificreports/ tianschanicus was almost straight compared to other two deer (Fig. 1). The cranial border of the scapula began at the cranial angle and extended cranioventrally, which was similar to the cranial border of H. i. argyropus. Slightly dorsal to the level of the scapular spine tuber, it ran caudoventrally toward the scapular neck to form the scapular notch. Then, it continued cranioventrally to form the supraglenoid tubercle. The caudal border, located along the line connecting the caudal angle to the caudal margin of the glenoid fossa, was oriented straight toward the scapular neck; moreover, it was thicker than the cranial border. A tuberosity produced by the attachment of the teres minor muscle was present on the scapular neck. A nutrient foramen was located on the caudo-proximal margin of the tuberosity. The subscapular fossa had a shallow concavity along the entire medial surface (costal surface). A round and low ridge was present at the caudal area of the fossa, forming the caudal scapular border (Fig. 1). The serrate surface extended from the cranial end to the caudal end of the medial surface, had a wide space at the cranial end, and became narrow caudally. The cranial and caudal serrate surfaces bordered the subscapular fossa along its lower ridge (Fig. 1).
Although scapular cartilage was lost during skeletal specimen preparation, there was clear evidence of its existence. Cartilage ossification is found in C. pygargus and H. i. argyropus aged> 4.5 years. The ossified cartilage is attached to the caudo-dorsal area of the cranial angle (Fig. 1). The cranial angle was round because the cranial border arose cranially, whereas the caudal angle was sharp and clear ( Fig. 1). There were no significant differences in the ventral angle among the ruminant species. However, a notch between the glenoid cavity and supraglenoid tubercle (NBSG) was absent from 12 C. p. jejuensis (12/18; 66.7%). The coracoid process was clearly observed in all species (Table 1, Fig. 2). The rate of glenoid notch appearance differed among species. The glenoid notch was commonly observed in H. i. argyropus (71%) and C. p. tianschanicus (78.3%), whereas it was absent from C. p. jejuensis (44.4%) ( Table 2, Fig. 2). No synovial fossa was present in the deer.
Osteometric data. The osteometric data are presented in Table 3. With the exception of Ld (longest dorsal length), the osteometric data significantly differed among species.The osteometric parameters include HS (height), DHA (diagonal height), Ld, SLC (shortest length of the scapular neck), GLP (longest length of articular process), LG (length of the glenoid cavity), BG (breadth of the glenoid cavity). The Ld was shorter in H. i. argyropus than in the other species, whereas it was similar between C. p. tianschanicus and C. p. jejuensis ( Table 3). The LTS (length of supraglenoid tubercle) was significantly shorter in C. p. jejuensis than in the other species. The SIs (scapula index) were 61.7%, 56.5%, and 56.1% in C. p. jejuensis, C. p. tianschanicus, and H. i. argyropus, respectively.
The result of discriminant functional analysis (DFA) showed that 97.2% of Korean wild deer were clustered accurately within each species (Fig. 4). Among 72 specimens, 2 C. p. jejuensis were misclassified as C. p. tianschanicus. About 91.6% specimens were assigned correctly in cross-validation procedure (p< 0.01) (supplementary data). The discriminant function 1 completely seperated H. i. argyropus and C. p. tianschanicus, but C. p. jejuensis overlapped with the other groups. On the other hand, the discriminant function 2 separated C. p. jejuensis from C. p. tianschanicus and H. i. argyropus groups. Among the predictors loadings on the discriminant function 1, HS, SLC, DHA, and BG were predictors which mainly contributed to seperating H. i. argyropus and C. p. tianschanicus, and DHA and BG contributed in it. The variables, HS, DHA, Ld, and LG, largely involved in discriminant function 2 which seperated C. p. jejuensis from the other species.
We compared the measurements between males and females (Table 3). In C. p. tianschanicus, the HS, Ld, SLC, GLP, LG, and BG values differed between sexes (HS, Ld, SLC, GLP, LG; p < 0.05, BG; p < 0.01). In H. i. argyropus, only Ld and SI differed between sexes (SI; p < 0.05, Ld; p < 0.01). No significant differences were observed between sexes in C. p. jejuensis.   21 findings of body size differences among these species. Body size can be influenced by various habitat conditions, and the variation in HS and Ld, which directly reflects overall scapular size, appears to be the consequence of adaptation to environmental conditions. The ratio of HS and Ld (SI) was higher in C. p. jejuensis compared to the other deer species ( p <0.01). This means that the scapula of C. p. jejuensis is relatively wider compared to other scapulae. Matsuo et al. 26 found that ungulates living in mountainous areas have wider scapulae than those living in open habitats. This finding correlates with the habitat preference of C. p. jejuensis, which is mainly found in oreums (small volcanic cones or hills that are formed by the accumulation of volcanic material around a vent) and Halla Mountain on Jeju Island 28 . Therefore, the distinct SI value of C. p. jejuensis seems to reflect its adaptation to the environmental conditions of Jeju Island.
The supra/infraspinous fossa is the attachment site of the supra/infraspinatus muscle and plays a significant role in determining muscle features 1,2,26 . Because the morphology of these fossae mainly depends on the position of scapular spine and the cranial border of the scapula, C. p. jejuensis which has a gentle "S" shape of cranial border of scapula possesses relatively larger supraspinous fossa compared to the other two species that has the straight cranial border of scapula. Conversely, the infraspinous fossa of C. p. tianschanicus is relatively larger than that of H. i. argyropus and C. p. jejuensis. The differences of morphology of supra/infraspinous fossae are associated with the supra/infraspinatus muscle features and the their locomotion capacity 26 . These morphological features may be induced by adaptations to their mountainous habitats and locomotion performance on sloped areas after considering the function of supra/infraspinatus muscle and the deer's habitat.
The scapular spines of Korean wild deer share several features with other ruminants. A weak scapular spine tuber, a thickened free margin of scapular spine were observed. The scapular spine arises from the dorsal border and curves caudally at the scapular spine tuber, which is similar to the morphologies in other ruminant species (except gorals) [8][9][10]13 . On the other hands, we found that the scapular spine slightly differed among species at the acromion level: the scapular spine was slightly inclined caudally in H. i. argyropus, whereas it was inclined cranially in C. p. jejuensis and C. p. tianschanicus (66.7% and 69.6%, respectively) ( Table 1). Hildebrand 10 found that scapular spine straightness could not be used to distinguish among species. However, previous studies [29][30][31] showed that mechanical stress on the skeleton causes skeletal development during prenatal and postnatal periods. Therefore, the different inclines of scapular spine can be affected by animal's habitual behaviors and adjusting to their habitats. The cranial acromion curve may result from the strong actions of supraspinatus and omotransverse muscles and the caudal curve could indicate that the caudal muscles attached to the scapula (i.e., thoracic region of the trapezius muscle, shoulder region of the deltoid muscle, and the infraspinatus muscle) are relatively larger in these taxa.
Although the development of scapular cartilage and its ossification with age in ruminant animals have been well-documented 7,8 , there is a lack of studies on the location and timing of this ossification. For instance, in horses, scapular cartilage appears after birth and fuses to the scapula after 3 years 8 . Wißdorf and Butendieck (1988) 13 identified the initial structure that becomes ossified in adult pudu deer shoulder cartilage; the structure also becomes attached to the cranial angle. In this study, bone fragments were observed in scapular cartilage mainly in deer aged over 4.5 years, suggesting that scapular cartilage ossifies with age in H. i. argyropus, C. p. jejuensis, and C. p. tianschanicus. The bone fragments in scapular cartilage exist at the cranial angle and near the caudal angle or dorsal border, with fragments near the cranial angle fusing with the scapula earlier than others. However, further research is required to determine the number of bone formation centers (ossification centers) in scapular cartilage, the timing of their appearance, and the timing of their fusion with the scapula.
The tuberosity associated with the attachment of the teres minor muscle, was clearly observed on caudal border of scapula in H. i. argyropus, C. p. jejuensis, and C. p. tianschanicus. Although this structure is evident in Bovidae animals, it is not frequently mentioned in books 7-9 , International Committee on Veterinary Gross Anomorphic Nomenclature 32,33 , or previous studies 10,13,34-36 . However, Nickel et al. (1986) 9 meontioned ridges on the caudal border close to the glenoid cavity in domestic ruminants, whereas Peters and Brink 37 observed this structure in two members of the Bovidae (springbok and grey rhebok).
In most animals, supraglenoid tubercle is connected to the glenoid cavity without a clear boundary although there is another ossification center that forms the tubercle and the cranial portion of glenoid cavity. However, in the present study, we noted a distinctive difference in scapulae of H. i. argyropus, and C. p. tianschanicus in terms of the presence of a clear boundary between the glenoid cavity and supraglenoid tubercle. The clear boundary formed the NBSG as shown in Table 2. Although these differences may not necessarily lead to differences in muscle attachment, they may be useful for species identification.
The glenoid cavity is the portion of the scapula that articulates with the humerus, and in ruminants, its margin is generally circular or elliptical, according to several studies 3,8,9,37 . While synovial fossa is present in some species like cattle, sheep, and goral, it is not found in others, such as pudu, fallow deer 13 , fossils of Rangifer tarandus, and Cervus elaphus 3 . In the current study, synovial fossa was not noted in scapulae of Korean wild deer.
Glenoid notch is a distinct structure which usually referred in horse osteology 8 . In ruminants, the notch is present depending on species such as goral 38 , buffalo 35 , and pudu deer 13 . We observed the glenoid notch in Korean wild deer. The rate of glenoid notch presence was lower in C. p. jejuensis than in H. i. argyropus and C. p. tianschanicus ( p. tianschanicus were similar (both 56%), whereas the SI of C. p. jejuensis was bigger than other two deer (62%). Thus, the scapula was lower triangular shape in C. p. jejuensis than in H. i. argyropus or C. p. tianschanicus. This is because the Ld did not differ between C. p. tianschanicus and C. p. jejuensis, whereas the HS differed between species. This difference may be related to the smaller body size and shorter leg length of C. p. jejuensis, compared with C. p. tianschanicus 21 . Although SIs are difficult to compare among species because these values were not consistently assessed in previous studies, Getty 8 remarked an SI of 60% in cattle; Hildebrand (1955) 10 found SIs of 62.5% in white-tailed deer, 69.0% in sheep, and 70.0% in goats; and Singh and Johari (1974) 36 found SIs of 45% in cattle and 56% in buffalo. The scapulae of H. i. argyropus and C. p. tianschanicus are similar to the scapulae of buffalo, whereas the scapulae of C. p. jejuensis are similar to the scapulae of white-tailed deer.
Sex differences were observed in the Ld and SI of H. i. argyropus and in the HS, Ld, SLC, GLP, LG, and BG of C. p. tianschanicus; however, no sex differences were observed in any measurements of C. p. jejuensis. The SI of the male H. i. argyropus is larger than the female because of a disparity in Ld scale. There is no sex difference in the SI of the C. p. tianschanicus. Notably, other measurments except for SI, LTS, and DHA of C. p. tianschanicus showed sexual dimorphism, but C. p. jejuensis did not. While deer exhibit sexual dimorphism in antlers and canines, there is no any study reporting sexual dimorphism in scapula. However, previous studies [18][19][20] have revealed scapular sexual dimorphism in some animal species, and there were different features between species. Therefore, the presence of sexual dimorphism can be used for interspecies distinction between C. p. jejuensis and C. p. tianschanicus.
In conclusion, we analyzed scapular sizes and shapes in H. i. argyropus, C. p. jejuensis, and C. p. tianschanicus, and found differences among the three deer. Morphologically, C. p. jejuensis demonstrated various differences compared to C. p. tianschanicus including the shape of the cranial border, the rate of glenoid notch presence and the presence of the NBSG. H. i. argyropus showed different acromion direction. In all osteometric data except LTS and SI, C. p. tianschanicus showed largest followed by C. p. jejuensis and H. i. argyropus. In SI, the Jeju roe deer was larger and the other two were similar. Sexual dimorphism was noted in SI and Ld of H. i. argyropus and 6 parameters (HS, Ld, SLC, GLP, LG, and BG) of C. p. tianschanicus, but not in C. p. jejuensis.These findings may provide an anatomical basis for the categorization of C. p. jejuensis and C. p. tianschanicus into subspecies.  Figure  LG, BG, and LTS. The data-set was standardised and used to establish classification of wild deer species. To verify the robustness of the classification, leave-one-out cross validation was used in R software.

Morphometric data collection.
For data citations of datasets uploaded to e.g. figshare, please use the SPSVERBc1 option in the bib entry to specify the platform and the link, as in the SPSVERBc2 example in the sample bibliography file.